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SUMMARY 


The origin of domestic dogs is poorly understood 
[1-15], with suggested evidence of dog-like features 
in fossils that predate the Last Glacial Maximum [6, 9, 
10, 14, 16] conflicting with genetic estimates of a 
more recent divergence between dogs and world- 
wide wolf populations [13, 15, 17-19]. Here, we 
present a draft genome sequence from a 35,000- 
year-old wolf from the Taimyr Peninsula in northern 
Siberia. We find that this individual belonged to a 
population that diverged from the common ancestor 
of present-day wolves and dogs very close in time to 
the appearance of the domestic dog lineage. We 
use the directly dated ancient wolf genome to recali- 
brate the molecular timescale of wolves and dogs 
and find that the mutation rate is substantially slower 
than assumed by most previous studies, suggesting 
that the ancestors of dogs were separated from pre- 
sent-day wolves before the Last Glacial Maximum. 
We also find evidence of introgression from the 
archaic Taimyr wolf lineage into present-day dog 
breeds from northeast Siberia and Greenland, 
contributing between 1.4% and 27.3% of their 
ancestry. This demonstrates that the ancestry of pre- 
sent-day dogs is derived from multiple regional wolf 
populations. 


RESULTS AND DISCUSSION 


The closest living relative of domestic dogs is the gray wolf, 
Canis lupus [1], but the number of domestication events, as 
well as their antiquity and geographical origin, is highly conten- 
tious [1-15]. While molecular estimates of the time of origin of 
the dog lineage are contingent on principally unknown mutation 
rates and generation times, the most recent genomic estimates 
of the divergence between wolves and dogs date to 11,000 to 
16, 000 years ago [13, 15, 17-19]. These estimates are in consid- 
erable discord with reported archaeological evidence of dog-like 
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canids from before the Last Glacial Maximum, which date as far 
back as 36,000 years before present (BP) [6, 9, 10, 14, 16]. 
Furthermore, a recent study showed that gray wolves from as 
disparate locations as China, Israel, and Croatia were symmetri- 
cally related to modern-day dogs [15]. This observation suggests 
that dogs were domesticated prior to the diversification of 
present-day gray wolf populations or that the wild ancestors of 
dogs are now extinct. The latter scenario would be consistent 
with an earlier finding of a morphologically distinct wolf popula- 
tion adapted to megafaunal prey in Late Pleistocene Beringia 
[20], as well as mitochondrial DNA evidence for a Holocene 
replacement of European gray wolves [21]. One hypothesis 
could thus be that the wild ancestors of dogs were a genetically 
distinct wolf population that inhabited the Late Pleistocene 
steppe-tundra biome and that this population was subsequently 
replaced [18], possibly by a northward postglacial expansion of 
smaller-bodied wolves that gave rise to modern-day wolf diver- 
sity. To test this hypothesis, we sequenced a draft genome of a 
Late Pleistocene wolf from northern Siberia. 

Shotgun sequencing was performed on a canid rib originally 
collected during an expedition to the Russian Taimyr Peninsula 
in 2010, from here on referred to as Taimyr 1. This rib was iden- 
tified as coming from a wolf through sequencing of a short region 
of the mitochondrial 16S rRNA gene and was directly dated using 
accelerator mass spectrometry (AMS) radiocarbon dating to 
30,920 + 380 '4C years BP, equivalent to ~34,900 calendar 
years BP (35,000 years ago) (for details, see the Supplemental 
Experimental Procedures). Genomic libraries were constructed 
with and without treatment with the uracil-specific excision re- 
agent (USER) enzyme mix [22]. This enzymatic treatment was 
employed in order to remove postmortem-derived uracil resi- 
dues that introduce errors into ancient genome datasets. We 
sequenced these libraries on the Illumina HiSeq platform to a 
total average 1-fold sequencing depth (Figure S1), the majority 
of which is derived from the USER-treated libraries. Using the 
retained postmortem damage patterns at methylated CpG sites, 
we observed nucleotide misincorporation patterns in these se- 
quences expected from DNA tens of thousands of years old (Fig- 
ure S1) [23]. We assembled the mitochondrial genome to an 
average 182-fold sequencing depth and reconstructed a phylog- 
eny relating the Taimyr 1 individual’s mitochondrial DNA lineage 
to those of modern-day dogs and wolves, as well as previously 
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Figure 1. Mitochondrial DNA Phylogeny 

The outgroups (one coyote and two Chinese wolf sequences) were excluded. 
Nodes with posterior probability >0.9 are marked with asterisks, and stars 
indicate ancient samples morphologically described as “dog like.” See also 
Figure S1 and Table S1. 


published mitochondrial DNA sequences from ancient canids 
[14]. The mtDNA of Taimyr 1 forms a distinct lineage separated 
from the dog lineages (Figure 1), in agreement with the generally 
basal grouping of all ancient wolves that has been observed in 
previous studies [14, 20]. The sequencing depth of chromosome 
X is 50.4% of that expected of an autosome of similar size, 
demonstrating that the Taimyr individual was male. 

A previous study used seven present-day canid genomes to 
show that all studied wolves share a common origin that ex- 
cludes the ancestors of domestic dogs [15]. The divergence be- 
tween the wolf and dog lineages was estimated to date back to 
11,000-16,000 years ago, but this estimate relied on assump- 
tions on the mutation rate, which has not been directly estimated 


for domestic dogs. A prediction of this model is that the ~35,000- 
year-old Taimyr individual would have lived long before the 
divergence of dogs and gray wolves, and thus would be sym- 
metrically related to both populations. In a principal component 
analysis (Figure S2), the Taimyr individual is approximately inter- 
mediate to gray wolves and dogs, but this could be expected just 
from the extra genetic drift that has occurred in those lineages 
since the death of the Taimyr individual. To formally test the hy- 
pothesis that Taimyr 1 lived prior to the divergence between 
dogs and gray wolves, we used D statistics [24, 25], expecting 
symmetry (D = 0) under the null hypothesis corresponding to 
the population history (Andean fox, (Taimyr 1, (dog, wolf))). This 
test, as well as all other population genetic analyses reported 
here, used the USER-treated portion of the sequence data in or- 
der to minimize the effect of postmortem degradation. It should 
be noted that we also found consistent results for an even more 
conservative set of analyses restricted to transversion polymor- 
phisms. Overall, we find that the data are consistent with the Tai- 
myr wolf being symmetrically related to wolves and dogs for all 
pairs of present-day canids (|Z| < 3). Conversely, tests with the 
Taimyr individual being assigned to a clade with one of the pre- 
sent-day canids to the exclusion of another, corresponding to a 
history (Andean fox, (canid,, (Taimyr 1, canids))), were all rejected 
(Z > 5). 

To investigate the phylogenetic position of the Taimyr wolf 
further, we fitted a population model to the genomes of the Tai- 
myr wolf and modern canids. We find that gene flow between 
present-day dogs and wolves after their initial divergence is 
required to explain the data (Supplemental Experimental Proce- 
dures). Once this gene flow is included in the model, the Taimyr 
wolf can be successfully fitted as belonging to the wolf lineage, 
the dog lineage, or the lineage ancestral to both wolves and 
dogs (Figure 2; Supplemental Experimental Procedures). In 
these three models, the Taimyr wolf is consistently placed very 
close to the split between the dog and wolf lineages, as 
measured in units of genetic drift (an Fs; of 0 to 0.007). Thus, 
our data are consistent with a trifurcation of the dog, wolf, and 
Taimyr lineages, indicating that they all diverged at about the 
same time. This appears to be inconsistent with the hypothesis 
that dogs and wolves diverged only 11,000—16,000 years ago, 
since under this model it might be expected that the Taimyr 
wolf would be confidently placed on the lineage ancestral to 
the dog-wolf split, due to the substantial amount of genetic drift 
that most likely would have occurred since the death of the 
Taimyr wolf (85,000 years ago) until the split between dogs and 
wolves some 20,000 years later. 

To estimate the divergence time between the Taimyr individual 
and present-day canids directly, and accounting for changes 
in effective population size, we estimated the probability 
F(Aderivea|Bheterozygous) [25] of an individual A (such as the Taimyr 
individual) sharing a derived allele discovered as a heterozygote 
in a diploid present-day individual B. At sites where the Chinese 
wolf is heterozygous, the Taimyr individual carries the derived 
allele at 30.8% + 0.1% of the sites, compared to ~32% for the 
Croatian wolf, dingo, and boxer (Figure 3; Table S2). We used 
this as a summary statistic to estimate the divergence time of 
the Taimyr lineage given a model of population history of the 
Chinese wolf inferred using the pairwise sequential Markovian 
coalescent method [26]. We find that calibration using the most 
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Figure 2. The 35,000 Year-Old Taimyr Individual Belonged to a Population that Was Genetically Close to the Ancestor of Present-Day Gray 


Wolves and Domestic Dogs 


(A) D statistics testing the phylogenetic relationship between the Taimyr wolf and present-day canids. The Andean fox is used as an outgroup in all tests. Error bars 


correspond to one block jackknife SE on each side. 


(B) A model of population history (admixture graph) fitted to the data. The Andean fox, which was used as an outgroup, is not shown. Although the graph shown 
here indicates minor gene flow from domestic dogs into regional wolf populations, we note that an alternative graph with gene flow from the wolf lineages into the 
domestic dog lineages provides equally good fit. In all fitted graphs, the branch length between the divergence of the Taimyr wolf lineage and the wolf-dog split is 
very short. Branch lengths are measured as the allele frequency variance along each branch (genetic drift) rescaled to be proportional to Fst. 


See also Figures S2 and S3. 


commonly assumed mutation rate of 1 x 10~® per generation and 
a 3-year gray wolf generation time [5, 15, 18, 19, 27] would imply 
that the Taimyr wolf diverged from the Chinese wolf 10,000- 
14,000 years ago (Figure 3), which is incompatible with its cali- 
brated direct radiocarbon date of ~35,000 years BP. Instead, 
the mutation rate must be substantially slower in order to be 
compatible with the age of the Taimyr individual, and we find 
that the Taimyr divergence can be accommodated by a mutation 
rate of 0.4 x 10~® per generation (Figure 3). However, it should be 
noted that this assumes that the Taimyr wolf is directly ancestral 
to the Chinese gray wolf. If there was structure between the an- 
cestors of the Chinese wolf and the Taimyr wolf, the mutation 
rate would have to be even slower, and as such a rate of 0.4 x 
10-® per generation is conservative. We emphasize that this mu- 
tation rate is for non-CpG sites, since SNPs in CpG dinucleotide 
context were excluded from the variants called in the present- 
day genome. Alternatively, our results could indicate that the gen- 
eration time is longer than 3 years, or Some combination of slower 
mutation rate and a longer generation time. Regardless, this 
direct evidence suggests a longer timescale of wolf-dog popula- 
tion history and thus implies that the 11,000—16,000 years ago 
wolf-dog divergence inferred in a previous study [15] should be 
recalibrated to ~27,000—40,000 years ago. 

To examine shared ancestry between the ancient Taimyr wolf 
and a larger set of modern-day dog populations, we used data 
from 48 dog breeds genotyped at ~170,000 SNPs [28] and 
computed D statistics to assess whether each breed shared 
more alleles with the Taimyr wolf than a set of 15 modern-day 
gray wolves (Figure 4; Table S3). We found clear evidence of a 
closer relationship between the Taimyr wolf and the Siberian 
Husky (Z = 4.3, p = 0.000009), Greenland Sledge Dogs (Z = 3.6, 
p=0.00016), and, to alesser extent, Chinese Shar-Pei and Finnish 
Spitz (p < 0.05) compared to other dog breeds (Supplemental 


Experimental Procedures). To estimate the proportion of ancestry 
derived from the Taimyr wolf lineage in the Greenland Sledge 
Dogs, we fitted an admixture graph using the Andean fox, pre- 
sent-day gray wolves, and German Shepherds. We find that the 
best-fitting graph posits 3.5% Taimyr-derived ancestry in the 
Greenland Sledge Dogs but that an ancestry proportion ranging 
from 1.4% and 27.3% is consistent with the data (Supplemental 
Experimental Procedures). These results can be explained either 
by a very early presence of dogs in northern Eurasia or by the ge- 
netic legacy of the Taimyr individual being preserved in northern 
wolf populations until the arrival of dogs at high latitudes. 

Extending our population history modeling to SNP array geno- 
types from a large set of gray wolf populations [1 1, 29], we further 
find that the majority of the ancestry of North American wolves 
also diverged from other wolves later than the Taimyr wolf line- 
age (Supplemental Experimental Procedures; Figure S4). This 
suggests that all extant gray wolf populations share a relatively 
recent origin, most likely sometime after the divergence of the 
Taimyr wolf lineage but prior to the inundation of the Bering 
Land Bridge and subsequent isolation of Eurasian and North 
American wolves. 

In conclusion, our results provide direct evidence for a longer 
timescale for the divergence of the dog and wolf lineages than 
previously assumed, and thus suggest that dogs may have orig- 
inated much earlier than commonly accepted. Such an early 
divergence is consistent with several paleontological reports of 
dog-like canids up to 36,000 years old [6, 9, 10, 14, 16], as well 
as the evidence that domesticated dogs most likely accompa- 
nied early colonizers into the Americas [30]. However, the initial 
divergence between the ancestors of dogs and gray wolves 
would not necessarily have had to coincide with domestication 
in the sense of selective breeding, since this human-mediated 
process could have occurred later or over an extended period 
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Figure 3. Recalibrating the Lupine Mutation Rate using the Directly 
Dated Taimyr Genome 

We estimated the probability F(A|B) for the Taimyr individual carrying the 
derived allele at positions where the Chinese wolf is heterozygous using 
empirical data. This statistic is informative on the coalescent time passed in the 
ancestry of the Chinese wolf since its divergence from the Taimyr wolf’s 
lineage. Since the Taimyr wolf must have separated at least 35,000 years ago, 
the age of the specimen obtained from a direct radiocarbon date, its diver- 
gence from the Chinese wolf can be used to calibrate the lupine mutation rate, 
in the sense that we can infer the maximum mutation rate that is consistent 
with the proximity of the Taimyr genome to the present-day Chinese wolf. To 
achieve this, we built calibration curves for F(A|B) given a model of Chinese 
wolf effective population size history (inferred using pairwise sequentially 
Markovian coalescent [PSMC] analysis) and different mutation rates. Mutation 
rates per generation are shown in the legend. A generation time of 3 years is 
assumed. See also Table S2. 


of time. On the other hand, a scenario of a much more recent 
timing of domestication would require that the majority of pre- 
sent-day dog ancestry originates from an extinct or presently un- 
sampled wolf population. Regardless, we find that the ancestry 
of present-day dog breeds descends from more than a single 
domestication event, since high-latitude dog breeds such as 
the Siberian Husky and Greenland Sledge Dogs can trace part 
of their ancestry to the now-extinct Taimyr wolf lineage. This 


introgression could have provided early dogs in high latitudes 
with phenotypic variation beneficial for adaptation to a new chal- 
lenging environment. 
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Figure 4. Introgression from a Population Related to the Ancient Taimyr Wolf into Northeast Siberian and New World Arctic Dog Breeds 
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comprised of 83 Boxer, Gordon Setter, Doberman Pinscher, and Newfoundland individuals. The color scale represents Z scores where the obtained D statistic 
values were normalized by SEs estimated using a block jackknife. Positive statistics suggest an excess affinity to the Taimyr wolf over the modern-day gray 
wolves for dog breed X. Negative statistics suggest an excess affinity to present-day gray wolves over the ancient Taimyr wolf for dog breed X. We tested 48 
different dog breeds with information from ~66,000 SNPs. The affinity of the Siberian Husky and Greenland Sledge Dogs to Taimyr is replicated also in statistics of 
the form D(Andean fox, Taimyr; present-day gray wolves, dog breed X). See also Figure S4 and Tables S3 and S4. 
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